Slave-Boson Three-Band Model with O-O Hopping for High-Tc Superconductors 
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Slave boson mean-field approximation is carried out analytically for weakly doped Cu02 con- 
duction planes, characterized by Cu-O charge transfer energy Apd, Cu-O hopping to, O-O hopping 
t' and repulsion Ud between holes on Cu site taken as infinite. At zero doping 5, finite negative 
t',|t'| < to/2, expands the range of stability of the covalent, conducting state on the expense of 
the insulating state which, however, remains stable at larger Ap^- For sufficiently large Apd the 
renormalized charge transfer energy saturates at A\t'\ instead of decreasing to zero, as at t' = case. 
In contrast to t' , finite 6 suppresses the insulating state nearly symmetrically with respect to the 
sign of 5. The regime with charge transfer energy renormalized close to 4|t'| fits remarkably well the 
ARPES spectra of Bi2212 and LSCO, and, in the latter case, explains the observed strong doping 
dependence of the Cu-O hopping. 



ARPES measurements of the electron spectra in the 
high-Tc superconductors have shown that the char- 
acteristic band energy scales fall into the range between 
0.1 and 1 eV. On the other hand, the high-energy spec- 
troscopies indicate considerably larger characteristic en- 
ergies, of the order of several eV This implies that 
the renormalization of the electron spectrum in low en- 
ergy range measured by ARPES is strong. Recent results 
on La2-a;Sra;Cu04 (LSCO) have shown in addition 
that the shape of the Fermi surface varies strongly with 
doping, indicating again strong renormalization of band 
parameters. Further argument in favor of the strong 
renormalization comes from the approximate electron- 
hole symmetry of the phase diagram of the high-Tc super- 
conductors If the couplings were weak, band renor- 
malization small, the approximate electron-hole symme- 
try would occur with respect to the (optimal) doping 5^, 
necessary to bring the Fermi level to the logarithmic van 
Hove singularity. This would be accompanied by the ac- 
tivation of the Umklapp scattering, i.e. tendency towards 
Mott localization. Thus far such tendency has not been 
observed and the approximate symmetry occurs with re- 
spect to (5 = rather than with respect to 5c- This rules 
out weak coupling and raises the question whether in- 
termediate 1^ to strong [0 coupling can account for such 
behavior. 

It is usually assumed that the largest coupling is the 
local hole-hole interaction Ud on the Cu-site of the con- 
ducting Cu02 planes. The interaction problem is often 
treated within the Emery model ||^ which describes the 
three-site structure of the Cu02 plane by copper and oxy- 
gen site energies ea and Sp respectively, the Cu-O hopping 
to and the 0-0 hopping tg, the single-particle parame- 
ters in addition to Urf. The strong renormalization of the 
band structure characterized by e^, Sp, Iq and tp is thus 
expected only when the number of holes on the Cu-site 
is close to unity, because only then large Ud is effective. 
In the opposite case, when the holes tend to reside on 



oxygen sites, the renormalization should decrease. Bear- 
ing this in mind, the limit of large is discussed here, 
taking Ud larger than Apd^ Sp — Ed, to, t'o. The limit 
of infinite is, as usual, treated within mean-field slave 
boson approximation pM, neglecting the AF and SC ef- 
fects. The latter occur [dT on the energy scales of the 0.01 
eV, an order of magnitude below the band energy scales 
of interest here. 

Previously, similar slave boson calculations were car- 
ried out analytically for the to, tg 1^ (CuO) analog of 
the 2d (CUO2) case [^, and, for the latter, numerically 
for some particular choices of the bare parameters P-pd]| . 
The 2d three-band structure is considerably richer than 
its Id two-band counterpart, therefore, similar holds for 
the corresponding mean-field slave boson theory. In par- 
ticular, the 2d case is sensitive to the sign of tg. Two 
bands may anticross, as in the Id case, but in 2d case, 
band-touching is also possible [ pd|jT^ ]. Above all, in con- 
trast to Id model, the parameter tp removes in 2d case 
Fermi energy of the half-filled band from the van Hove 
singularity, i.e. introduces a finite doping Sc required to 
bring it back ]T^ . The question whether the electron- 
hole symmetry occurs for S — or 6c should therefore be 
discussed at finite t'o. 

The slave boson procedure amounts to searching the 
ground state band energy Eg at a given number of holes 
1+S, subjected to the restriction that the number of holes 
Ud on the Cu-site is smaller than 1, as required by large 
Ud- In the mean-field approximation Eg is the energy 
of the free holes with the renormalized band parameters 
A,t, t', chosen to fulfill the requirement (rid) < 1, which 
becomes Rl 



{nd)+f/tl^l, (nd) 



dEp 
dA ■ 



(1) 



The minimization of Eq with respect to the remaining 
parameter t, gives the second slave boson equation 
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(ub) = 2 



t(Apd - A) 



{riB) = 



dEo 
dt ' 



(2) 



denoting renormalized bond charge on the Cu-0 bond. 

The solution of the coupled integro-differential equa- 
tions (|l]) and (|^) gives A and t as functions of Ap^;, 
t' and 5 (reniembering that if, = t'). The present slave 
boson calculation is carried out for arbitrary Ap^ and 
to, assuming 2\t'\ < to and \d\ « 1. The assumption 
2\t'\ < to is in the spirit of the Emery model but 
the opposite limit will nevertheless be briefly commented 
upon. On the other hand, the results obtained here for 
I (5 1 << 1 do not cover all the physically achieved doping, 
but rather indicate the tendencies associated with finite 
S. 

It is possible to transform jl^ the integro-differential 
equations (|^) and (||) to the coupled algebraic equations 
for all interesting regimes of parameters by calculating 
Eo{A,t,f , S) analytically and then (nd) and (?t._b) from 
Eqs. (]^) and (||). Leaving the details of the calculation 
for the extended publication jl^ , the final results will be 
discussed here. 

These results are best understood starting from the 
t' = situation. Sp, which does not renormalize, is cho- 
sen for the energy origin, Sp = 0. Only one band, s^{k), 
out of three is occupied for S < 1, 
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/i = sin^ 



16t2/i), 



2 ky 



(3) 



For Apd close below A 
expression t. 



setting A{t' — 0) = Aq). Using e° (k) calculate Eo, 
leads to the solution of the t'=0 slave boson equations 
(|l|) and (|^) for arbitrary Ap^ and to (taken positive with- 
out loss of generality). The behavior of t for (5 = 
is fairly simple: t/to = 0, starting from large Ap^ to 
A'Zi = 4.7tn, where the angular point in t/to occurs. 

t is given by a Landau-like 

A^^ = l.lStoy^A^^- Ao, which describes 

Brinkman-Rice (BR) phase transition from the insulat- 
ing i = to the conducting phase. The overall behavior 
can be obtained by expanding Eq.(H) in terms of t^/Ao 
including the terms of the order of t close below A^^ . 

The corresponding behavior of Ao(Aprf) is more impor- 
tant for further discussion. Aq has a single, maximum 
(angular at 5 = 0) at Apd = A^^, where it reaches the 
value Aq = Ap^/2. For Apd well below A^^ the regime 
t > \Ao\ is reached, (n^) decreases, according to Eq.(|^) 
and so do both renormalisation of to to t and of Ap^ to 
Aq. Before approaching asymptotically Ap^ for the large 
negative Apd /to, Aq crosses zero at Apd /to ~ 1.5: 



standard covalent situation Ap^ 
extreme Aq ~ covalent limit. 



to renormalizes to the 
For A > A-, A(Apd) 



decreases fast towards the asymptotic regime valid for 
large Ap^, when A = 4.7t§/Apd tends to zero (similarly 
to the site energy of the t-J model JT^). 



Turning next to finite t', in the limit where 2\t'\ < to, 
it is to be noted that under such condition (more pre- 
cisely for A\t'\ < A^y2 = 2Ato), 4|i'| intersects Ao(Apd) 
at A' > A-. When \t'\ « to. A' = 4.7i2/|t'|. For 



Apd < Ap^, A\t'\ is smaller than Aq and/or t obtained 
from t'=0 solution, and a perturbative treatment of t' is 
possible. For Ap^ > A^^, 4|i'| > Aq and full treatment 
of t' is required. 

The perturbative corrections due to t' in the regime 
t > Aq at Apd well below 4.7^0, are not of particular in- 
terest because only of the quantitative nature. The atten- 
tion will be therefore focused here on the small t regime 
at Apd close below or above 4.7to- For this purpose a 
perturbative calculation can be carried out in terms of 
small t, using (with Sp = 0) 



ei(k) = -A - At' 



o{t'). 



A/i - St'h 

A2 - 16t'2/2 

=sin2^sin2^. 
■' 2 2 



(4) 



This expansion, valid at (5 ~ for < t < A — A\t'\ and 
arbitrary sign of t' , reduces at t'=0 to the spectrum used 
in obtaining the BR transition, i.e it is appropriate for 
following the evolution with t' of the small t slave boson 
solution. 

In the first step Eq.(U) is used in the range A.lto < 
Apd < Ap^, where it can be expanded further in terms 
of < |t'|/A < 1. The following results are obtained re- 
taining the terms linear in t' and quadratic and quartric 
in t. 

The corresponding analytical calculation of the po- 
sition of the Fermi level sp has been carried out ear- 
lier | p2[ , by systematical linearization with respect to 
t'. The shift of the ep from the van Hove singularity 
at Ao{evH — £f) = 2tjjXF turned out to be 



, 1 87r2(,5, -,5) 
a;F m — = 77 — 

Xp 1 + OTT^Oc 



(5) 



where 5c = —?>2t'/TT'^Ao stands for the critical doping 
required to bring the Fermi level back to the van Hove 
singularity. Experimentally, this doping is positive i.e. 
signi' — — signA. The linearization procedure can be 
generalized to the calculation of Eo{A, t, t', 5), and then, 
straightforwardly, to the solution of the slave boson equa- 
tions (|l|) and (||). Here, for brevity, small terms in Eo, 
linear in xp, discussed in Ref. p^] ,will be omitted in the 
description of this solution. 

The BR transition is shifted to the new position A^^ — 
A.7to + 3.5Sc, with — 13.6t'/to, but it is not suppressed 
by t'. This important result can be traced back to the 
outstanding feature of the expansion of Eq. (^) in terms 
of t'/A, generating the term t'^t'/A'^. The above term 
vanishes with t, making the theory all along convergent 
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with respect to i'/^j exemplified by Eq.(||), where xp 
is finite and small independently on t. 

Acr, the value of A in the angular point, is equal 
to the Ap^/2, i.e. for 5c > it is sHghtly increased 
with respect to the t' = value, and A remains above 
Aq ~ A{t' = 0) in the whole range of convergence of this 
approach,A^^ < Ap<j < A^^. 

These results can be extended to finite S > 0. In 
contrast to t' , i.e. to 5c, finite S suppresses the BR 
transition, turning it into crossover, in spite of Eq.(^, 
which suggests analogous roles of dc and S. In the 
vicinity of A^^, the angular behavior of A is essentially 
modified by finite S. For example, the value of A at 



pel 



4.7io + 3.55c is 



^ - A = 0.8(1 + 0.2(5c)5^/^ 
and correspondingly 

t^/tl = 0.63(1 - 0.26c)S^^^, 



(6) 



(7) 



up to the leading orders in S and 6c- 

Finally we turn to the range Apd > A^^, where the 
stability of the t = solution is to be discussed using 
Eq.(^ non-expanded in terms of t', but dropping out 
O(t^) term. The main contribution to the cohesive en- 
ergy Eq is weakly dependent on £p (i.e. Sc) and strongly 
dependent on lower cut-off, which is A = A — 4|t'| ac- 
cording to Eq.(^, i.e. Sc- (Note in this respect that the 
e F dependency of Eq was not very important even in the 
vicinity of the BR transition, Eqs.(^and (^.) Expanding 
then Eq.(^ around k = (tt, tt) point, one finds immedi- 
ately that Eq is logarithmic in t'^/t'A, which, through 
slave boson equations, leads for Ap^ :s> A'^ at (5 = to 



A = A-4U'| w 2\t'l 



(8) 



where c is a numerical factor. Eq.(H) shows that the 
asymptotic behavior of A which at t' = vanishes as 
Aq ~ ip/Apc;, undergoes an essential modification to the 
exponential behavior, with the saturation at A — |4i'|. 
Although the derivation of the result (H) does not suggest 
that this result is independent of the dimension, the sim- 
ilar behavior was obtained |^ at 5 = in Id. Moreover, 
it is easy to see that finite S > leads to t ^ ^/S as be- 
fore. However, the convergence of this result is reduced 
to t < A, which is itself (exponentially) small according 
to Eq.(H). The regime i « A is quickly reached with 
doping and it is clear that, from the slave boson point of 
view, the limit t > A is physically achievable. 

The corresponding fits of the ARPES data for the 
Fermi surfaces of LSCO at different doping d are shown 
in Fig.l. These fits also require signA = — signi' and de- 
termine two parameters e.g. A/t' and t/t' as functions of 
6. While A w 4\t'\ is almost independent on doping, t in- 
creases quickly with doping. For the underdoped regime 



(small doping) the behavior of t can be approximated 
by \/5, corresponding to the heavy-fermion limit, as sug- 
gested by here-presented mean-field slave boson theory. 
However, for larger doping, appreciable values of t are 
reached, leading the system to the essentially covalent 
behavior. Finally, for overdoped regime, t and A seem 
to saturate, requiring further theoretical investigation, in 
particular because the Id theory ||^ suggests otherwise. 
Here, the calculus becomes quite complicated since the 
perturbative treatment of t is no longer possible, and the 
exact hole spectra are required. Fig.|^ shows the con- 
ducting band fit of Bi2Sr2CaCu208 (Bi2212), which fixes 
the absolute value of t in addition to the relative values 
obtained from the Fermi surface fits, t turns out to be 
of the order of O.leV, almost an order of magnitude less 
than to, again in agreement with the present slave boson 
analysis. 
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FIG. 1. (a-e) Experimentally measured Fermi surfaces for 
LSCO |l[ (green dots) and the three band model fits (solid), 
(f) fitting parameters in units of t' with A small 

Equally good fits (c.f. Fig.||b) were however obtained 
earlier Q assuming |A| < 4|t'| rather than |A| < 4|i'|. 
Since it was shown here that the latter is the solution of 
the (5 « slave boson theory when 2\t'\ < to, the for- 
mer can only occur at 5 « 0, if 2\t'\ > tg. While the 
behavior of A given by Eq.(j^) can be visualized, assum- 
ing for the moment i = 0, as the avoided crossing of 
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the dispersionless (t=0) copper level and the pure oxy- 
gen band with the energy — 4|i'| at k = {tt,tt) point, the 
limit |A| < 4|i'| means that the Cu level has entered the 
oxygen band. The holes on Cu site are then necessar- 
ily transferred to the oxygen sites, i.e. (n^) < 1. For 
n = I + S, with (5 ~ 0, Eq.(|^) then requires finite t, i.e 
the solution t=0 is unstable for 2\t'\ > ig- (For finite t, 
copper band undergoes anticrossing (sign t' = - sign A), 
rather than touching (sign t' = sign A) with the oxy- 
gen band. With (n^) < f, the renormalization 
becomes small, analogously to the case of small and neg- 
ative Apd at t'=0. t < to and A < Apd are not expected 
to depend much on doping S. The appreciable exper- 
imental dependence of t on small 6 [0 is therefore an 
indication in favor of the strongly renormalized 2\t'\ < to 
solution. 




FIG. 2. Band-fits for the conducting band of Bi2212 of |||, 
(a) A small, (b) A small 

Let us finally point out that 2\t'\ < to solution, strongly 
dependent on doping, is approximately symmetric with 
respect to the electron {S < 0) or hole doping {6 > 0). 
As the latter was discussed above, leading in particu- 
lar to t/to ~ VS for sufficiently large Ap^ > 4.7to, 
the case S < will be briefly mentioned here. Since 
(ud) < n = 1 + 5, Eq.(||) means that t > toV—^, i-e-, 
roughly speaking, t behaves as ^/\5\, irrespectively of the 
sign of ^. t ~ line of the slave boson theory, obtained 
for the half-filled S — Cu-based band, is the singular 
line of the phase diagram, not only at t'=0, but also for 
t' 0, as shown here. In the properties depending on 
the details of the band structure close to the Fermi sur- 
face, the simplest of which is perhaps the Hall constant, 
the symmetry is broken when a finite t'is taken into ac- 
count. The corresponding asymmetry is the largest in 
the small t'/A limit of Eq.(||), when i5 > means doping 
towards and S < away from the vH singularity, whereas 
in the opposite limit of Eq. , the Fermi level at ^ = 
is so far from the vH singularity, that the variation of 
the topology of the Fermi surface with doping becomes 
small. Overall, the asymmetry is much smaller than in 
the weak coupling limit. 

In conclusion, the present work shows that the band 



structure of the high-Tc superconductors can be under- 
stood in terms of p-d bands in the Cu02 plane, renormal- 
ized by strong interaction on the Cu- site. According to 
the low energy ARPES data, the bare band parameters, 
measured by high energy spectroscopy, are renormalized 
by an order of magnitude. More profound understand- 
ing of the energy scales, of the order of O.I eV, i.e. the 
understanding of the approximate ground state on these 
energy scales, represents an appropriate prerequisite for 
consideration of the symmetry breaking, associated with 
AF and SC, on the energy scales of 0.01 eV. 
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